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ABSTRACT 



Context. Infrared Dark Clouds (IRDCs) harbor progenitors of high-mass stars. Little is known of the parental molecular clouds of the 
IRDCs. 

Aims. We demonstrate the feasibility of the near-infrared (NIR) dust extinction mapping in tracing the parental molecular clouds of 
IRDCs at the distances of £) » 2.5 - 8 kpc. 

Methods. We derive NIR extinction maps for 10 prominent IRDC complexes using a color-excess mapping technique and NIR data 
from the UKIDSS/Galactic Plane Survey. We compare the resulting maps to the '^CO emission line data, to the 8 pm dust opacity 
data, and to the millimeter dust emission data. We derive distances for the clouds by comparing the observed NIR source densities to 
the Besanfon stellar distribution model and compare them to the kinematic distance estimates. 

Results. The NIR extinction maps provide a view to the IRDC complexes over the dynamical range of Ay « 2 - 40 mag, in spatial 
resolution of ~ 30". The NIR extinction data correlate well with the '^^CO data and probe a similar gas component, but also extend to 
higher column densities. The NIR data reveal a wealth of extended structures surrounding the dense gas traced by the 8 pm shadowing 
features and sub-mm dust emission, showing that the clouds contain typically > 10 times more mass than traced by those tracers. The 
IRDC complexes of our sample contain relatively high amount of high-column density material, and their cumulative column density 
distributions resemble active nearby star-forming clouds like Orion rather than less active clouds like California. 
Conclusions. NIR dust extinction data provide a new powerful tool to probe the mass distribution of the parental molecular clouds 
of IRDCs up to the distances of Z) ~ 8 kpc. This encourages for deeper NIR observations of IRDCs, because the sensitivity and 
resolution of the data can be directly enhanced with dedicated observations. In addition to mass distributions, the NIR data provide 
relatively reliable distance estimates. 

Key words. ISM: clouds - evolution - ISM: structure - dust, extinction 



1. Introduction 

The initial conditions of molecular gas from which high-mass 
stars form remain poorly established. This is, for the most part, 
due to the violent nature of high-mass star formation which, once 
started, impacts the cloud structure on a very short timescale. 
Consequently, the sites where these initial conditions might still 
be prevailing can only scarcely be found from the nearby high- 
mass molecular regions such as the Orion, Cepheus, or Vela 
clouds where the star formation clearly has already taken place 
(and is still on-going). 

Infrared Dark Clouds (IRDCs) were originally identified as 
dark structures a gainst the brig ht mid-infrared background of 
the Galaxy (.Perault etal.lll996t lEgan et al.|[T998l) . Due to the 
details of this technique, the IRDCs do not form a homoge- 
nous sample of molecular clouds, but preferably high contrast 
(~ high-column density) objects at the distances of D < 8 
kpc are detected. The IRDCs are objects in which conditions 
suitable for high-mass star and star clu ster formation are be- 
lieved to be present or developing (e.g., iRathborne et al.ll2006t 



IVasvunina et al.l 120091 : IRathborne et al.1 l2010h . This picture is 
based on the fact that they contain massive clumps of cold 
gas (T ^ 15 K), up to the masses of young cluster- f ormiri g 
clumps (e.g., ICarev et al.1 119981: IRathborne et al] l2006l l2010l) . 
Furthermore, these clumps can contain substructures, i.e. cores, 
with sizes comparable to pre-stellar low-mass cores, but with 
significantly higher masses (and hence, densities). Some IRDCs 
have been observed to harbor candidate intermediate- or high- 



mas s protostars (e.g., iP illai et aL[ 
Beu ther et alJl2010l: lHen ning et aL 
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; IRathborne etaP 120101: 
I). This demonstrates the 
capability of IRDCs to form high-mass stars, although to what 
extent is still under debate (e.g., Kauffmann & Pillai 2010). 

The observational technique through which the IRDCs are 
identified (mid-infrared shadows against the Galactic emis- 
sion) is sensitive to the material in the column density range 
mUj) ^ 10-100 x 10^^ cm-Me.g- iButler & TanI 120091: 
iPeretto & Fulled 120091: IVasvunina et alJ l2009l) . Similar sensi- 
tivity is re ached by mapping thermal dust emission from 
IRDCs (e.g..'Rathborne et al."2006^ 'Beuther & Steinacke ijl2007l: 
iRathborne et al. 2010; Schuller et al. 2009). Mid-infrai-ed ex- 
tinction mapping using background sources can reach slightly 
better sensitivity, but the resolution remains worse due to low 
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density of background sources (108", cf. iRygl et alJ l2010h . 
Observations of CO emission lines can reach relatively low col- 
umn densities, thereby p robing the envelopes of mid-infrared 
IRDCs (e.g.. iHernandez & Tan 201 li), but bec ome optically 
thick, or freeze out (e.g.. iHernandezet al.ll20Tll) . at A^(H2) > 
5 - 10 X 10^' cm"^. As a result, information on the large- 
scale {I > a few pc) environment of IRDCs, and of the connec- 
tion of such component to the IRDCs, has remained virtually 
non-existent. Consequently, some basic questions regarding the 
parental clouds of IRDCs remain unanswered. To what degree 
are the IRDCs close to each other in fact components of coher- 
ent, large-scale filaments or other structures? Are the molecular 
clouds harboring IRDCs similar to the nearby high-mass clouds? 
How much material there is surrounding IRDCs and what is the 
role of that material for their evolution? 

Near-infrared (NIR) observations of stars that shine through 
molecular clouds provide an alternative technique to probe the 
low column density, cloud-scale structures surrounding IRDCs. 
The NIR dust extinction mapping technique is commonly used 
to estimate gas colum n densities towards local (r < 1 kpc) 
molecular clouds (e.g., 'Lada et al."1994'; 'Cambresv et al."2002( 
Kainulainen et al. 2006, 2009b; Lombardi et al. 2006, 2008a, 
20101: IRowles & Froebrichi 2009). Such techniques are based 



on measuring either the interstellar reddening towards sources 
that are behind the molecular cloud (color- excess met hods, e.g. 
iLada et alJ[T994l: iLombardi & AlvesllIOOll: lLombardill2009l) . or 
the surface density o f such sources (star-count methods, e.g. 
ICambresv et al.ll2002h . The benefit of such techniques is two- 
fold. First, the relation between NIR extinction (or redden- 
ing) and dust column density depends only relatively weakly 
on the optical properties of dust grains. Furthermore, it does 
not depend on dust emissivity index or on dust temperature. 
Second, the extinction mapping in NIR is sensitive over a con- 
siderably different column density range, i.e . ~ 10^' N < 



6 X 10^^ cm~^ (e.g.,O mbardi & Alves 20 0TI: iKainuIarnen et al.l 
l2007t ILombardi et"aiT T2008a: Roman-Zuniga et al.ll2010h . than 
the mid-infrared shadowing and dust continuum emission obser- 
vations. Compared to mid-infrared observations, it is also un- 
affected by local variations in the cloud's surface brightness, 
thereby being relatively unaffected by signatures of star forma- 
tion such as nebulae. These properties make near-infrared dust 
extinction mapping an important complement to mid-infrared 
extinction and mm dust emission data. 

In this paper, we examine the feasibility of using the NIR 
extinction mapping method to systematically probe the parental 
molecular clouds of IRDCs that are located at the distances 
of several kiloparsecs at the Galactic plane where the confu- 
sion from diffuse dust is significant. We derive extinction maps 
for a sample of 10 prominent IRDC complexes us ing archival 
NIR d ata from the UKIDSS/Galactic Plane Survey (Luc as et all 
l2008l hereafter GPS) and compare them to '^CO emission line 
data from the Boston University-FCRAO Galactic Ring Survey 
jjackson et al.,,20 06, hereafter GRS), mid-infr ared opacity data 
from t he recently published IRDC catalogue bv lPeretto & FuUen 
(I2009h . and to sub-millimeter dust continuum data from the 
ATLASGAL survey (Schulleret al. 2009). This combination of 
data will be used to examine the mass distribution surrounding 
the IRDCs and the fraction of dense cold gas in their parental 
molecular clouds. We will also use the NIR data to derive 
distances for our IRDC sample by determining the density of 
sources on the foreground to the IRDCs and comparing that to 
the predicted sourc e density from the Besan5on Galactic stellar 
distribution model dRobin et al.ll2003h . 



In Section ^we introduce the data used in the paper and the 
extinction mapping technique. In Section [3] we present the dis- 
tance determination method and apply it to the sample. Section 
^ presents our results, which are then discussed in ^ In 
Section |6]we give our conclusions. In addition. Appendix lAl 
gives the extinction maps for the sample, together with the '^CO, 
8 jum opacity, and sub-mm dust emission data. 



2. Methods and data 

2. 1 . Sample of IRDC complexes 

We selected the objects for this study by visually inspecting 
the Spitzer/GLIMPSE survey images and choosing from them 
a sample of spatially extended complexes with high contrast 
against the background. The sample selection was restricted, on 
the one hand, by the coverage of the UKIDSS/GPS that was used 
as NIR data in the study, and on the other hand, by the cov- 
erage of the GLIMPSE survey and the contrast of complexes 
in it. These restrictions resulted in selection of complexes be- 
tween / ^ - 40 deg and \b\ < 1 deg. One of the main objectives 
of our study is to test the practical limits of the NIR extinction 
mapping technique in mapping distant clouds. Therefore, we se- 
lected complexes located both relatively near (e.g., 35.49-00.31, 
Dg„ = 3.0 kpc) and far (e.g., 32.09+00.09, Dg,,s = 7.1 kpc) based 
on kinematic distance estimation from the GRS data. As a re- 
sult, the sample does not form a well-defined selection, although 
the chosen clouds are likely among the largest (most massive) 
IRDC complexes at the chosen spatial range. Nine out of ten 
sample clouds have counterparts in the 8 fj.m opacity catalogue, 
and eight in the BU/GRS survey. All complexes are included in 
the area surveyed by ATLASGAL. The chosen complexes are 
listed in Table [T] and designated accord ing to the notation of the 
GRS survey (Ro man-Duval et al.ll20Tol) . 

2.2. Near-infrared extinction mapping 

We employed the color-excess mapping technique nicer 
(Lombardi & Alves 2001) to derive dust extinction maps for the 
sample of IRDC complexes. The method was used in conjunc- 
tion with JH Ks band photometric data from the UKIDSS/GPS 
(iLucas et al.l 12008). In particular, data from the Data Release 7 
Plus were used. The data reach approximately the limiting mag- 
nitudes of K - 19 mag, H = 19.1 mag, and 7 - 20 mag. The 
UKIDSS project is defined in Lawrence et all (l2007h. UKIDSS 



uses the UKIRT Wide Field Camera (WFCAM iCasah et al 
2007). The photometric system is descri bed in |Hewett et al 
(2006), and the calibration is described in Hodgkin et al. 
The science archive is described in Hambl y et al. ( 2008) . 

In the basic implementation of the nicer extinction mapping 
technique the NIR colors of stars, shining through molecular 
clouds, are compared to the colors of stars in a nearby control 
field that is free from extinction. This comparison yields mea- 
surements of a NIR reddening towards the stars in the molec- 
ular cloud region, from which extinction can be estimated by 
adopting a wavelength dependency for the reddening, i.e. a red- 
dening law. The extinction values derived towards each source 
(Ay) are then used to compute a spatially smoothed map of 
dust extinction through the cloud (Ay). This basic approach is 
straightforward to apply for nearby molecular clouds (D < 500 
pc), in which case almost all sources that are detected are lo- 
cated behind the intervening dust cloud. However, in the case of 
IRDCs which are known to be at the distances of several kilo- 
parsecs, a significant fraction of the detected stars are located 
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between the cloud and the observer, i.e. on the foreground to 
the cloud. Such stars bias the extinction determination at low 
colu mn densities and totally hamper it at high column densities 
(e.g.. lLombardil2005h . Optimally, the foreground sources should 
all be eliminated from the sources used for extinction mapping. 
In practice, excluding them uniformly from everywhere is not 
possible due to the spread in the intrinsic colors of stars which 
mimics dust extinction. In the following, we introduce the prac- 
tical implementation of the nicer technique used in this paper 
For the descriptio n of th e basics of the method itself, we re- 
fer to iLombardi & AlvesI (12001) (see also Lombardi 2005, who 
presents an analysis of the effect of foreground stars for the tech- 
nique). 

We retrieved from the UKIDSS archive the JHKs data for 
the chosen complex es. For the N IR reddening-law, we adopted 
the coefficients from lCardelh et a l. (1989): 



TK = 0.600 X = 0.404 X Ty = 0.114 X Ty. 



(1) 



The extinction mapping procedure itself consisted of three steps 
which were required especially to deal with the large number 
of foreground sources. This procedure is illustrated as a flow- 
chart in Fig. [1] In the first step, a "dirty" extinction map was 
calculated in order to find a low-extinction region to be used as 
a control field for the next step of the mapping. The map was 
also used to determine high extinction regions from which the 
surface number density of foreground stars, needed in the next 
step of the map derivation, could be determined. 

In the second step, the contribution of foreground stars was 
statistically subtracted from the mean colors of the control field 
in order to calculate the mean color of stars on the background 
to the IRDC complex. In the statistical subtraction, the density 
of the foreground stars in JHK color-color space was first sub- 
tracted from the density of all control field stars to yield a fore- 
ground corrected source density in each color-color bin. Then, 
the source number in each color-color bin of the control field 
was reduced to correspond to this corrected value. The fore- 
ground corrected mean colors were then computed using this 
reduced set of control field stars. The subtraction is important, 
because the control field is not a truly non-extincted field, but 
suffers from extinction due to the diffuse, extended dust compo- 
nent along the line of sight in the Galactic plane. This diffuse ex- 
tinction at the plane v aries between Ay « 1 - 4 mag kpc"' (e.g., 
[Marshall et alJl2006l) . As a result, the stars farther away than 
an IRDC complex are expected to be clearly redder than those 
closer than it. In order to estimate the reddening caused by the 
IRDC complex itself, it is then necessary to estimate the color 
of the background population by eliminating the contribution of 
the foreground sources from the mean colors of the control field. 
The reddening of the sources in the control field is illustrated in 
Fig. 121 left panel, which shows the NIR color-color diagram of 
the sources in the control field chosen for 38.94-00.46. For bet- 
ter visibility, only sources with the photometric errors cr < 0.3 
mag are plotted. The figure clearly shows how a large fraction of 
sources is heavily reddened along the reddening line. 

The effect of removing the contribution of the foreground 
stars from the control field is demonstrated in Fig. |2] center 
panel. The figure shows the color-color diagram of all sources 
in the control field of the 38.94-00.46 complex with a greyscale. 
The greyscale represents the number density of stars in each 
color-color bin. The mean colors calculated for this distribution 
is shown with a green plus sign. The red contours in the fig- 
ure show the number density of the foreground stars identified 
from the 38.94-00.46 cloud region (identification procedure is 
explained in more detail later). The colors of foreground stars 
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Fig. 1. Flow diagram illustrating the implementation of the 
nicer technique in this paper (explained in Section 92.2b . The 
acronyms^ and bg refer to foreground and background, respec- 
tively. 



are clearly less reddened, as expected for stars that on average 
are closer to the observer and thus suffer less from diffuse ex- 
tinction. Using the surface number density of these foreground 
stars, we subtracted their statistical contribution to the mean col- 
ors. This resulted in estimates of the mean colors of the popu- 
lation on the background to the cloud complex, indicated in the 
figure with a green cross. This change in mean color was signifi- 
cant for all mapped clouds, and clearly required for a reasonable 
definition of the zero-point in the extinction mapping technique. 

The variability of the diffuse extinction component (and the 
uncertainty in determining the foreground population in the con- 
trol field) induces uncertainty to the zero-point of the extinction 
determination. We assessed the level of this uncertainty by per- 
forming for one cloud the zero-point determination using 10 dif- 
ferent control fields. The standard deviation of the derived zero- 
points was about Ay - 1 mag. We note that such assessment 
was not possible for all clouds (or for numerous control fields 
for one cloud), since finding non-extincted control fields close 
to the complexes is a problem in general. We consider the result 
of this experiment to be indicative of the zero-point uncertainty 
in our method. 

In the third and final step the extinction map was corrected 
for the contribution of foreground sources in the field. The iden- 
tification (and removal) of foreground sources is trivial in re- 
gions where the extinction is high compared to the scatter in the 
intrinsic colors of stars. The uncertainty of the extinction deter- 
mination for an individual source is approximately cr(Ay) a: 3 
mag, so the sources that are located in a region of intermediate 
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local mean of extinction (9 mag < Ay < 20 mag mag) can be 
trivially removed. The effect of foreground stars can be elimi- 
nated also at lower extinctions by estimating the surface number 
density of foreground objects from high Ay regions and sub- 
tracting the estimated contribution of the foreground population 
to the local mean colors. This procedure is explained in the fol- 
lowing. 

We estimated the surface number density of the foreground 
population by examining the sources in regions where the local 
mean extinction is higher than Ay > 6 mag (about 2-cr eiTor of 
an individual extinction measurement). Fig.|2]right panel illus- 
trates this procedure by showing the frequency distribution of 
individual extinction measurements, A^, towards such regions 
of the 38.94-00.46 complex where the local mean extinction is 
above Ay > 6 mag. Clearly, the foreground stars form a sepa- 
rate 'bump' in the distribution slightly below Ay = mag. The 
surface number density of the foreground sources was estimated 
by fitting a Gaussian to this bump and integrating the number of 
stars within it. 

In principle, it would be possible to use the determined fore- 
ground source density to statistically subtract the foreground 
stars from the observed field down to Ay =0 mag. In practice, 
however, doing this would impose an assumption that all extinc- 
tion features are caused by a dust component at one well-defined 
distance. In order to avoid this assumption, we only performed 
the elimination of the foreground stars down to Ay regime in 
which the foreground star density could be determined, i.e. down 
to such Ay values in which the foreground bump could be iden- 
tified from a frequency distribution such as shown in Fig.|2]right 
panel. In practice, this requirement resulted in elimination of 
foreground stars above Ay > 4 mag. Thus, the extinction maps 
above this limit are, to first order, free from foreground star con- 
tamination, but suffer from it below it. However, we note that 
below this limit the foreground star contamination actually be- 
comes clearly lower, because the relative amount of the back- 
ground sources is high. 

The spatial resolution of the extinction data is set by the 
width of the Gaussian smoothing function that is used in smooth- 
ing the pencil-beam extinction measurements towards stars (Ay) 
onto a map grid (Ay). In practice, the choice for this is guided 
by the surface density of the background sources, as all map pix- 
els towards which extinction is to be calculated should have at 
least one actual measurement (source) inside. For the depth of 
the UKIDSS data, we concluded that the spatial resolution of 
30" represents a good balance between higher resolution and 
lower noise at intermediate column densities. This translates to 
the physical resolution of 0.4 pc at 2.5 kpc distance. Typically, 
the pixels with Ay = 20 mag have a few sources inside the beam 
area (depends on the Galactic coordinates). 



2.3. CO data 

We used the publicly available '^CO emission line data from the 
GRS as another tracer of gas in the mapped IRDC complexes. 
These data provide a tool to independently verify the physical 
connection between clouds located close to each other in projec- 
tion. The GRS survey provides '^CO data for the galactic lon- 
gitudes / = [18°,55?7], thus including 8 of the 10 complexes 
of our sample. The spectral resolution of the data is 0.2 km s 
the typical rms noise 0.4 K, and the spatial resolution is 46". 
From these data, Rathborne et al. (2009) identified structures in 
position-position-velocity space, resulting to a catalogue of 829 
molecular clouds containing 6124 clumps. We use the designa- 



tions of the GRS catalog as hsted in lRoman-Duval et al ] (l2010h 
to refer to the clouds of our sample (see Table[T]l- 

We constructed position-velocity diagrams for each region 
that was mapped in NIR extinction and used them to disentangle 
the gas component primarily responsible for the extinction fea- 
tures (see Fig.|3]for an example). This component was then cho- 
sen from the spectral line cubes and integrated to yield the total 
emission from the cloud complex. The maps of integrated '^CO 
intensity are shown for each complex in Appendix lAl (see also 
Fig. 13). In almost all chosen complexes, the dust extinction fea- 
tures clearly match the morphology of the CO data at the chosen 
velocity interval, thereby indicating that the extinction features 
indeed are caused by the gas component in the chosen veloc- 
ity interval. The correlation between NIR extinction and CO is 
further discussed in ^4.11 



2.4. Mid-infrared opacity data 

We compared the NIR extinction data also to the mid-infrared 
dust opacity maps derived from 8 jum shadowing feat ures in 
the im ages of the Spitzer/GLIMPSE survey by Peretto & FuUeii 
( l2009h . While that technique reaches higher extinctions (up to 
Ay < 100 mag) in high resolution (~ 1") compared to NIR 
mapping, it is considerably less sensitive to low extinction val- 
ues (Ay 1 - 10 mag). Therefore, the two techniques comple- 
ment each other well, with the NIR data probing lower density 
envelopes and the 8 fim data probing higher density cores. 

For each IRDC complex, we retrieved the 8 fj.m opac- 
ity mapfl of all indivi dual IRDCs identified in the region by 
Pere tto & Fulleii (|2009). Typically, there were several tens of 
these IRDCs within the mapped regions. From the maps of in- 
dividual IRDCs, we then constructed 8 jum opacity maps for the 
entire complexes. These maps are shown alongside the NIR ex- 
tinction maps and the '^CO integrated emission maps in Fig. [3] 
and in Appendix|A] In general, the correspondence between the 
NIR and MIR extinction data is rather scarce, due to both very 
different spatial resolution and sensitivity. The data are com- 
pared more closely in ^4.11 



2.5. 870 fim dust emission data 

Thermal dust emission from the Galactic plane (|/| < 60°, 
\b\ < 1.5°) has been mapped at 870 fim by the Atacama 
Pathfinder Experime nt (APEX) within the ATLASGAL survey 
(ISchuUer et al.ll2009l) . The survey provides emission data in res- 
olution of 19.2" with typical sensitivity of 0.25 Jy/beam (5cr). 
This sensitivity limit translates to A^(H2) ~ 6.1 X 10^' cm"^, 
or Ay 6.1 mag, using typical values of T - 20 K (e.g., 
Eganet al.1 19*9 8) and /fgvo = 1.85 g cm ( Ossenkopf & Hennina 
1994t) . The thermal dust emission observations will not "sat- 
urate" at high column densities, unlike NIR and MIR extinc- 
tion techniques, and thereby they provide perfect complimen- 
tary data to trace the densest parts of the IRDC complexes. The 
data from ATLASGAL for each complex are shown in Figs. |3] 
and AppendixlAlto allow comparison with the other data. In 14. 21 
we use these data to estimate the mass of dense gas in the com- 
plexes. 



3. Distance determination using NIR source counts 

Nine out of ten complexes in our sample are covered by the 
GRS survey and kinematic distances have been derived for 
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Fig. 2. Left: Color-color diagram of the sources in the control field chosen for the 38.94-00.46 complex with photometric error 
<T < 0.03 mag. Extinction of Ay = 5 mag is indicated with an arrow. Center: The same for all sources in the control field, shown 
with a grey-scale (and black contours) that indicates the density of points. The mean colors of these sources is marked with a green 
plus sign. The red contours show the density of the sources that are on the foreground to the 38.94-00.46 complex. The contribution 
of these sources to the mean color is statistically subtracted from the mean, resulting to the mean color of the background sources 
that is shown with a green cross. Right: Frequency distribution of the extinction measurements towards individual sources (Ay) in 
a region where Ay > 6 mag. The small bump at slightly negative extinction values are the foreground stars that are removed during 
the last step of the extinction mapping procedure (see Section ^2.2i . The surface density of the foreground sources identified in this 
way is also used in estimating the distances to the cloud complexes in Section ^ 



them bv lRathborne et all (l2006h andlRoman-Duval et al.l (l2009h . 
Using a different approach.'Mars hall et al.l (12009 ') used NIR data 
from the 2MASS survey together with the Besangon Galactic 
stellar distribution model to derive distances towards a large 
number of IRDCs, including six complexes of our sample. We 
list the distances from both sources in Table [1] using nota- 
tion Dpr s for kinematic distance measurements and Dgen for 
[Marsh all et al. | (|2009) distances. Similarly with the technique 
adopted by Marshall et alj (l2009h . the photometric data from 
the UKIDSS/GPS employed in this paper can be used to de- 
termine cloud distances. With the help of the derived extinc- 
tion maps, it is straightforward to measure the surface density 
of stars that are located between the cloud and the observer, i.e. 
on the foreground to the cloud. The measured surface density can 
then be compared to what is predicted for the Galactic coordi- 
nates of the cloud by the stellar distribution model of the Galaxy, 
whic h yields an estimate for the cloud distance (e.g.. lAlves et alj 
Il998i) . This technique have been demonstrated to yield dis- 
tance estimates that are in agreement with esti mates made with 
other techniques for nearby molecular clouds (Lada et al.l l2009l; 
[Lombardi et al. 201 1). We used this approach to estimate the dis- 
tance for the clouds in the sample. In the following, the distance 
determination procedure is explained in detail. 

For each mapped complex, we estimated the surface density 
of the foreground stars by counting the sources in K band in re- 
gions where the derived extinction was between 5-8 mag < 
Ay < 20 mag, depending slightly on the Galactic coordinates of 
the cloud (because the error of the extinction measurement de- 
pends on the coordinates). In these regions, the foreground stars 
are easily distinguished from background stars due to their low 
extinction values compared to the local mean extinction (illus- 
trated in Fig. 12] and explained in l2.2l i. For the distance determi- 
nation, it is preferable to count the foreground stars down to as 
faint magnitude as possible in order to keep the uncertainty due 
to counting statistics (Poisson error) as low as possible. Yet, the 
sources should only be counted down to the completeness limit 
of the photometric data, Kutn- Since the completeness limit of the 



UKIDSS/GPS data changes as a function of the Galactic coor- 
dinates, the limiting magnitude above which the source counts 
were performed changed from cloud to cloud. Typically, the 
sources brighter than Ku^n < 16.5 mag were included. Typically, 
the regions from which the sources were counted were ~ 50 
arcmin' (~ 800 map pixels) and contained ~ 1000 foreground 
sources. 

The observed foreground source densities were then com- 
pared to the source counts predicted by the B esangon Galactic 
stellar distribution modefl ( Robin et al.ll2003l) . The model was 
used to retrieve differential source counts in the distance bins of 
0.2 kpc, down to the previously chosen limiting magnitude /Tum, 
at the positions of the mapped clouds. For the purpose of dis- 
tance determination, the most crucial parameter in the stellar dis- 
tribution model is the amount of the diffuse, inter-cloud extinc- 
tion from which all sources, including the foreground sources, 
suffer. This diffuse extinction generally amounts to ~l-4 mag 
kpc ' and can vary rapid ly as a functio n of Galactic position 
(e.g.. [Marshall et alj|2006h . In particular. [Marshall et al.1 (l2006h 
used the 2MASS NIR data to derive radial extinction profiles in 
the Galactic plane in spatial intervals of 15' and radial bins of 
~ 200 - 500 pc. We used these data to estimate the level of dif- 
fuse extinction towards the IRDC complexes. This was done by 
choosing 2-4 positions surrounding the complex (depending on 
the morphology of the complex) and computing the mean ex- 
tinction profiles for those lines of sight. As an example. Fig. [5] 
shows profiles for two clouds, one showing a very linear dif- 
fuse extinction component (38.94-00.46) and another showing 
a more complicated diffuse component (34.24+00.14). The fig- 
ure shows two individual profiles towards these clouds with plus 
signs and dotted lines, and the mean of these profiles is shown 
with a red solid line. 

The differential source counts retrieved from the Besangon 
model were then modified (reddened) using the mean diffuse ex- 
tinction profiles shown in Fig. [5] These reddened count data were 
used to compute the final predicted source densities as a func- 



Available on-line at |http://model.obs-besancon.fr/| . 
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Fig. 5. Top row: Radial dust extinction profiles constructed for 
two complexes. The plus signs and the dotted lines show the ex- 
tinction profiles adopted from Marshall et al. (2006). The solid 
red line shows the mean profile. Bottom row: Predicted source 
count functions (source count above th e limiting magnitude 
/Tiim), based on the Besangon model and iMarshall et alJ (l2006t) 
reddening data, for the same two clouds. The source count func- 
tions corresponding to the mean extinction profile are shown 
with the solid red curves. The dashed black curves show the 
same function using the individual extinction profiles. The hor- 
izontal lines show the observed foreground source density and 
the vertical lines point out the corresponding distances. The dot- 
ted lines show the Poisson errors related to the observed source 
density (and the resulting distances). 



tion of distance. Then finally, the observed foreground source 
densities were compared to this function to yield the distance 
estimate. Fig.|5]illustrates the predicted cumulative source count 
functions as a function of distance. The source count function 
corresponding to the mean extinction profile is drawn with a red 
solid line. The observed foreground source density is indicated 
with a solid horizontal line, and the associated Poisson counting 
error that results from the finite region available for the source 
density determination is shown with dotted lines. This statistical 
error is typically » 100 pc. The resulting distance is shown with 
solid vertical lines. We refer in this paper to the distances de- 
rived from the foreground source densities with Dfg and list the 
distances for all clouds in Table [T] 

Figure |5] also shows the predicted source count functions 
computed using two individual extinction curves available for 
the cloud (dashed curves). This illustrates the uncertainty in the 
distance determination due to the uncertainty in the radial ex- 
tinction profile. For complexes with almost linear extinction pro- 
files, the scatter in the resulting distances is typically 5-10 %. For 
complexes showing more complex profiles it can amount up to 
~ 15 %. This clearly dominates the error of the distance mea- 
surement over the statistical uncertainty. We conclude that the 
uncertainty of our distance determination is about 15 %. 



4. Results 

4. 1 . IRDC complexes as viewed by NIR extinction 

We used the technique described in ^2.21 to derive dust ex- 
tinction maps for the 10 cloud complexes listed in Table [T] 
Figure[3]shows, as an example, the map derived for the complex 
38.94 - 00.46, together with the '^CO integrated line emission, 
the 8 yum dust opacity, and 870 yum dust emission maps. We give 
in AppendixlAlthe same figures for all complexes of the sample. 
In addition. Fig. |4] shows a smaller-scale zoom-in to the central 
region of the complex 35.49-00.31 to demonstrate the relatively 
high resolution of the NIR extinction data. 

In general, the NIR extinction maps reveal extended (tens-of- 
parsec size), highly structured cloud complexes in the mapped 
regions. In the chosen resolution (30", corresponding to 0.44 
pc at 3 kpc distance), the maps probe extinction values up to 
Ay < 30 - 40 mag. The maps commonly reveal how the IRDCs 
(identified from 8 fim data) in the regions are connected by 
lower column density (Ay ~ 5 mag) structures. While these 
low-column density structures already strongly suggest a physi- 
cal connection between the IRDCs in the regions, the connection 
is confirmed by the '^'CO data showing the structures to be con- 
nected in the position-position-velocity space. 

Generally, the NIR extinction features resemble closely the 
morphology of the '^CO emission in the chosen velocity inter- 
vals. To illustrate the correspondence. Fig. |6] shows a pixel-to- 
pixel comparison of the NIR and '^CO data for one of the clos- 
est clouds (38.94-00.46, D - 2.7 kpc) and for the farthest cloud 
(32.09-1-00.09, D = 8.0 kpc). The comparison shows a clear cor- 
relation between the quantities. In most cases, as for example in 
the case of 38.94-00.46 shown in Fig.|6] the CO emission satu- 
rates around Ay ^4-6 mag, a feature commonly observed in 
CO-to-Ay ratios of molecular clouds and associated to optical 
dept h effects and/or depletio n of CO at higher column densities 
(e.g.. lCombes 1991; Pineda et aLlliOOl . The figure also shows 
the ratio for these quantities measured typically at low column 
densities in local m olecular clouds (i n particular, for the Perseus 
molecular cloud, Pineda et an i2008h . At low column densities, 
the observed data roughly agree with this relation. 

The NIR extinction data and 8 fim opacity data have a factor 
of ~ 15 difference in spatial resolution and ~ 5 in sensitivity, 
which makes comparing them directly with each other difficult. 
Qualitatively, the morphology seen in the NIR extinction (and 
CO) maps is usually recognizable also in the 8 jjm cloud de- 
tections. This is clearly shown in, e.g.. Fig. |3] for the complex 
38.94-00.46 and in Fig. [Ajfor 11.10-00.10. This coiTelation 
is no surprise, because the complexes for the study were cho- 
sen as high contrast features from the Spitzer/GLIMPSE images 
(as described in 92. 11 1. Again qualitatively, the morphologies 
of smaller-scale structures, i.e., individual IRDCs, vary tremen- 
dously between the NIR and 8 yum data. While in many cases 
no IRDCs are detected at the positions of NIR extinction peaks 
(and vice versa), in other cases the correlation is very good. As 
an example. Fig. |4] shows a zoom-in to the central region of the 
35.49-00.31 complex. The 8 jum, NIR and CO detections in this 
region correlate very well. 

It is worth noting that unlike 8 yum data the NIR extinc- 
tion mapping is relatively unaffected by the signatures of star- 
forming activity in the clouds. An example of this is shown in 
Fig. [3] in which a Spitzer/GLIMPSE/MIPSGAL 3-color image 
of 38.94-00.46 is showrQ. The cloud region contains several neb- 
ulae that are likely excited by the formation of intermediate- or 

^ Available through [http://www. alienearths . org/glimpse/ 1 . 
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A, [mag] A, [mag] 

Fig. 6. Pixel-to-pixel comparison of NIR extinction and '^CO in- 
tegrated line emission for the 38.94-00.46 (Dfg - 2.7 kpc) and 
32.09-h00.09 (Z)fg = 8.0 kpc) complexes. The dotted line shows 
a typical ratio of these quantities for local molecular clouds 
JPineda et al...2008.) . 

high-mass stars. Clearly, these nebulae affect strongly the MIR 
opacity map, practically suppressing the detection of the MIR 
shadowing features ("holes" in the 8 yum opacity map). However, 
the nebulae are not as strong sources in the NIR and thereby do 
not strongly affect the NIR extinction map. 

4.2. Masses of the complexes from different tracers 

The NIR extinction data can be used to estimate the total gas 
column density along the line of sight by first converting the 
NIR extinction to visual extinction (Eq. [TJ, and then converting 
the the visual extinction to the total hydrogen column density 
using the measured ratio (i Bohhn et al. 1978) 

2N(}i2) + NiU) = 1 .9 X 1 O^" X Av [cm^^ mag" ' ] . (2) 

To calculate the total masses from the column densities, we 
summed up all the extinction values in the clouds above Ay - 
2.5 mag. This relatively high threshold was chosen to better dis- 
entangle the cloud complexes from the surrounding, possibly 
physically unrelated extinction components. We used the dis- 
tances for the clouds derived from the foreground stellar density 
measurements (see Sec. [3]), Dfg, that are listed in Table[r|in cal- 
culating the masses. The resulting masses are listed in Table [T] 
and they are generally between 1 and 10 x lO'* Mq. For compar- 
ison, the total masses of the nearby Taurus and Orion A molec- 
ular clouds above Ay = 1 mag are ~ 10'* M p and ~ 10^ M©, 
respectively (e.g.. iKainulainen et al.ll2009bl) . However, a large 
fraction of the total mass of the cloud is likely in regions below 
our selected threshold of Ay > 2.5 mag. This fraction is gener- 
ally ~ 4 - 80 %, and for example ~ 43 % in Taurus an d Orion 
A (e.s.. I Kainulainen et al.ll2009bt iLombardi et al.ll2010l) . Thus, 
most complexes included in this study have masses on the order 
of ~ 10^ Mq, equaling roughly the mass of the Orion A molec- 
ular cloud. We note that the largest uncertainty in the derived 
masses originates from the distance uncertainty. For example, a 
distance uncertainty of 500 pc for a cloud at 4 kpc distance re- 
sults in ~ 30 % uncertainty in total mass. The uncertainty of the 
total masses induced by the uncertainty of the extinction zero- 
point determination (Ay w 1 mag, see Section 32.2b is on the 
order of 15 %, resulting from the fact that typical mean extinc- 
tions above the mass calculation threshold are < Ay >^ 6-7 
mag. 

For comparison, we also list in Table [T] the total masses of 
the clouds calculated from the GRS CO line emission data by 



(iRoman-Duval et ani2010h . corrected to correspond to the dis- 
tance that was used in calculating the masses from NIR data. 
In the Roman-Duval et al. (2010) work, the cloud masses were 
calculated using '^CO line emission. The excitation tempera- 
ture was estimated from the '^CO line assuming that it is op- 
tically thick. This temperature was adopted for calculation of 
'^CO optical depths and thereby '^CO column densities. Finally, 
the H2 column densities were estimated from the '^CO col- 
umn densities using the re lations »('^CO)/»('^CO) = 45 and 
«('^C O)/«(H7) = 8 X 10-5 (i Langer & Pen^ll990HBlake et al.l 
[T987I) . The masses calculated from CO emission are relatively 
close to the masses derived from NIR extinction, although they 
seem to be systematically larger by a factor of ~ 1.5-2. This 
difference likely has its origin i n the different techniques u sed to 
identify clouds from data sets. 'Roman- Duval et all dloTO ) em- 
ployes thresholding in position-position - veloc ity space using the 
CLUMPFIND algorithm by Willian^ ([1991 . while we identi- 
fied clouds with a simple column density threshold in the plane 
of the sky. The difference in the total masses may also be in- 
duced by the fact that the zero-point determination of the extinc- 
tion mapping is not very robust for clouds at large distances (see 
Section 32.2l for discussion). 

We also calculated the total masses of the mapped regions 
from the 8 yum opacity data. For the transformation between 
opacity and column density, we adopted the reddening curve be- 
tween 8 fim and NIR as given by Indebetouw et al. (2005), i.e. 
T&fim = 0.4 3 X tk- With the redd ening curve between NIR and 
visual from lCardelli et all (Il989l) . this yields the conversion fac- 
tor Tg^m = 0.048 X Ay. Then, Eq. |2] was again used to trans- 
form Ay to hydrogen column densities. The resulting masses 
are listed in Table [T] The masses from 8 yum emission are typ- 
ically on the order of ~ lO'* M©, i.e. about an order of magnitude 
smaller than the NIR masses. This is quite expected, because 
the 8 fim data probe (only) significantly higher column densities 
than NIR data, thereby not detecting the bulk of the cloud mass. 
Obviously, the 8 fim absorption is spatially less extended than 
NIR or CO. This comparison of the total masses is given solely 
to demonstrate the extent of the diffuse component in the clouds 
not detectable via 8 |/m data. 

The masses of the dense gas component in the complexes 
were also calculated from the 870 fim dust emission data. In this 
calculation, we assumed for simplicity a constant temperature 
of r = 20 K for all complexes. Then, assuming optically thin 
emis sion, the colurnn den sities can be calculated from the equa- 
tion (ISchuUer et al.i2009l) 



' By{Td)Q.Ky^miH ' 

where R is the gas-to-dust-ratio, R - 100, By is the Planck func- 
tion at temperature T^, ^ is the beam size, Ky - \ .85 g cm"^ is 
the mass absorption coefficient at 870 //m, and ju = 2.8 is the 
mean molecular weight. The resulting masses are given in Table 
[1] The masses traced by the dust emission are generally on the 
order of ~ 10 - 20 % of the masses traced by NIR and CO data, 
and thereby close to the 8 //m masses. 

5. Discussion 

The NIR extinction mapping technique employed in this paper 
provides a view to molecular clouds harboring potential progen- 
itors of massive stars over a column density range Ay 2 - 40 
mag, with a spatial resolution of ~ 30". This range greatly com- 
plements other commonly used column density tracers, i.e., CO 
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Table 1. Cloud complexes and their properties. 



Complex 




D ^ 




MnirLIO'^ Mo] 






MsiO^m/MfiiR [%] 


03.80-01.00 


2.5 


- 


- 


1.7 


- 


- 


12 


11.10-00.10 


4.1 


4.74 


3.6 


4.5 


- 


11 


15 


18.54-00.16 


5.9 


4.33 


4.1 


4.7 


1.1 


17 


22 


24.94-00.16 


3.5 




3.4^ 


3.6 


2.7 


10 


10 


32.09 + 00.09 


8.0 


7.45 


7.r 


12 


1.5 


29 


26 


34.24 + 00.14 


5.9 




3.8 


13 


2.7 


46 


51 


34.39-00.71 


0.7 




1 


0.9 


0.3 





3 


35.49-00.31 


3.2 


3.63 


3.0 


5.6 


0.7 


8 


7 


36.74-00.16 


5.1 


3.88 


3.6 


5.4 


1.7 


7 


14 


38.94 - 00.46 


2.7 


4.14 


10.5 


2.0 


1.6 


50 


14 



" Distance derived from the foreground source densities in this paper (see Section[3]l. 

* Distance derived from modeling of NIR source color distribution bv ' Marshall et al.' ("2009*). 

Unless otherwise stated, the distance is a kinematic distance from the GRS survey (Roman- Duval et aDl2009h . 

As estimated in lRoman-Duval et alj j2010l) . corrected to correspond to the distance Dfg listed in column 2. 
" iRathbomeetal] ( 120061) 



line emission, mid-infrared dust extinction, and thermal dust 
emission. While CO obviously can provide kinematic informa- 
tion not reachable via NIR data, the dynamical range of NIR data 
extends to clearly higher column densities and is free from de- 
pletion/optical depth effects. This allows the NIR data to probe, 
in a uniformly calibrated manner, the regime from diffuse en- 
velopes to the denser clumps in the clouds. Compared to other 
dust tracers, the NIR data is more sensitive to low column den- 
sities. It can also provide a reasonable overlap with the dynami- 
cal ranges of dust emission and MIR extinction, allowing cross- 
calibration. Perhaps most importantly, the NIR extinction map- 
ping can provide temperature independent data on the mass dis- 
tributions of IRDC complexes in spatial resolution that matches 
that of the Herschel satellite that will provide high sensitivity 
dust emission data of the Galactic plane (e.g., iMoIinari et alj 
[2010). 

We independently derived the distances for the complexes 
by comparing the observed foreground stellar densities to the 
predictions from the Besangon stellar distribution model. We es- 
timated the uncertainty of our approach to be about 15 %, re- 
sulting mostly from the variable diffuse extinction component 
in the Galactic plane. There is a reasonable agreement between 
the kinematic distance estimates and earlier distance estimates 
based on modeling stellar distributions using 2MASS data (Z)gis 
and Dgen in Table [T] respectively), considering the uncertainties 
of each method. With only 10 clouds in our sample, we did not 
quantify the differences in the distances derived with different 
methods in more detail. Compared to the kinematic distance es- 
timates, the insensitivity of our approach to the near/far distance 
ambiguity is an asset. For example, in the case of the 38.94- 
00.46 complex it seems clea r that the kinematic distanc e deter- 
mination algorithm used by iRoman-Duval et a (120091) has in- 
correctly chosen the far distance for the cloud. As a downside, 
the foreground stellar densities can only be used to determine 
distances to clouds within ~ 8 - 10 kpc distance. In summary, 
the distance measurements confirmed the NIR extinction map- 
ping technique to be feasible for clouds as far as ~ 8 kpc away, 
thus allowing a vast sample of IRDCs to be mapped with the 
method. 

Compared to mapping local (D < 500 pc) molecular clouds, 
applying the NIR dust extinction mapping technique to clouds at 
the distances of a few kiloparsecs has two main obstacles. First, 
the large fraction of foreground stars towards the clouds biases 
the determination of extinction at Ay < 4 mag. As discussed in 



Section 32.21 below this column density the contribution from 
such stars biases the measurement towards lower Ay values. 
Second, the distribution of diffuse dust at the Galactic plane 
can have rapid spatial variations making determining the red- 
dening caused solely by the cloud in question uncertain. These 
two short-comings have the main effect that the zero-point of the 
Ay mapping is clearly less well-defined as in mapping nearby 
clouds. This can have a significant effect on the estimates such 
as the total masses of the clouds, since most of the molecular 
cloud mass is always at low column densities. 

The gas component traced by the NIR extinction matches 
closely the morphology of the CO line emission, with the dif- 
ference that it extends to clearly higher column densities (i.e.. 
Ay ^ 30 - 40 mag). Consequently, the masses of the com- 
plexes derived from NIR data are comparable to the ones de- 
rived from CO emission. Comparison of these masses with the 
masses derived from 8 jum dust opacity and 870 /urn dust emis- 
sion shows that on average these higher column density tracers 
reveal ~ 10 - 20 % of the masses of the clouds inside which they 
are embedded (with a large scatter). However, as speculated in 
Section l4r2l it is likely that our total mass estimate from NIR ex- 
tinction still significantly underestimates the total masses of the 
complexes due to the relatively high threshold for mass deter- 
mination, and thereby the masses revealed by 8 //m dust opacity 
and 870 fim dust emission are even clearly smaller, likely < 10 
% of the total masses. 

Figure |7] shows the cumulative mass functions calculated 
from the NIR extinction data. The functions indicate what frac- 
tion of the cloud mass is at column densities higher than the 
reference value (abscissa). For this plot, all the maps were 
smoothed to common resolution of 0.86 pc, reflecting 30" reso- 
lution at the 5.9 kpc distance (except 32.09+00.09 which is at 8.0 
kpc distance). For comparison, the figure shows the mass func- 
tions of the nearby Orion A and California mo l ecular c louds de- 
rived from similar data by iKainulainen et al.l (l2009bl) . Clearly, 
the IRDC complexes of our sample contain relatively large 
amount of high density material. Their mass functions seem to 
be markedly closer to that of the active star-forming Orion A 
cloud than the more quiescent California cloud. In particular, 
none of the IRDCs shows as steep mass function as California. 
We note that the mass functions presented here are based on 
data that have ~ 8 times lower spatial resolution than data of 
Kai nulainen et al.l ( l2009bl) . and thereby the relations may not be 
directly comparable. However, we find them clearly indicative. 
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Fig. 7. Cumulative mass functions for the complexes, indicat- 
ing the fractional mass above a chosen column density (or vi- 
sual extinction. Ay) threshold. The dotted red and dashed blue 
lines show the relation for Orion A (star-forming) and California 
(more quiescen t) clouds, respectively, both of which are about 
10^ Mq clouds dKainulainen et al.ll2009bi) . 

implying that the IRDC complexes show relatively high frac- 
tion of high-column density material comparable to very active 
nearby molecular clouds. 

While recent studies based on the molecular line data of 
the GRS have generated a census of molecular clouds in the 
Galactic plane (e.g. , Jackson et al. 2006; Rathborne et al. 2009; 
iRoman-Duval et al] f2010l) . the dust extinction maps can provide 
a mass tracer for these clouds extending over a clearly wider dy- 
namical range than the CO data. Most importantly, the dynami- 
cal range will overlap with that of higher column density tracers 
(especially thermal dust emission), allowing combining the in- 
formation from the two tracers. Such an approach can provide 
column density data extending from low column density enve- 
lope material to gravitationally dominated objects in the clouds. 
Thereby, the data sets together can be used to reveal the de- 
tailed mass distributions surrounding the progenitors of high- 
mass stars, and we will continue towards deriving such data in 
forthcoming work. 

6. Conclusions 

In this paper, we examined the feasibility of NIR dust extinction 
mapping technique in tracing low-to-intermediate column den- 
sity structures surrounding prospective birthplaces of high-mass 
stars, i.e., IRDCs. We used the data from the UKIDSS/Galactic 
Plane Survey to derive dust extinction through 10 cloud com- 
plexes at the distances between ~ 2.5 - 8 kpc, harboring alto- 
gether hundreds of IRDCs. We compared the derived NIR ex- 
tinction maps to the '^CO molecular line data from the Boston 
University-FCRAO Galactic Ring Survey, 8 yum dust opac- 
ity data from a recently published Spitzer IRDC catalogue by 
[Peretto & Fuller (2009), and to 870 ^mi dust emission data from 
the ATLASGAL survey dSchuller et al.ll2009l) . The conclusions 
of our work are as follows. 

1 . We demonstrated that the NIR dust extinction mapping tech- 
nique provides a powerful tool to trace the mass distribu- 
tion of typical IRDC complexes. Using only moderately deep 
UKIDSS/GPS data (K < 19 mag), extinction maps covering 



the dynamical range of Ay 2 - 40 mag can be constructed 
in the resolution of FWHM x 30". Such data provide highly 
complimentary, temperature independent mass distribution 
data to be combined, e.g., with data soon available from the 
Herschel satellite mission. Similarly, the results encourage 
for dedicated, deep NIR observations of IRDC complexes, 
because the spatial resolution of the extinction data at the 
high source density of the Galactic plane can be greatly im- 
proved with deeper observations. 

2. The NIR extinction maps reveal the low column density 
(Ay < 10 mag) structures surrounding IRDCs in great de- 
tail. The IRDCs are dense clumps in these parental (giant) 
molecular clouds that contain typically ~ 5 - 10 % of the 
total masses of the clouds. 

3. Despite their low star-forming activity the IRDCs of our 
sample contain relatively large fraction of high column den- 
sity material, as judged by their cumulative mass distribu- 
tions. Their mass distributions resemble more closely that 
of active star-forming clouds (like Orion) than less active 
clouds (like California). 

4. We derived distances for the complexes in our sample by 
comparing the observed NIR foreground source counts to 
the Besangon Galactic stellar distribution model. We esti- 
mate the accuracy of the distance measurements to be about 
15 %. The distance measurements confirm the feasibility of 
the extinction mapping method up to distances of ~ 8 kpc. 
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Fig. 3. Infrared Dark Cloud complex 38.94-00.46 at the distance of Dfg = 2.7 kpc (estimated from foreground source density, see 
Section ® . Top row: The first panel shows NIR dust extinction map, derived using JHKs data from the UKIDSS Galactic Plane 
Survey. The resolution of the map is 30" and contours starting from Ay = 3 mag in steps of 4 mag are shown. The second panel 
shows the integrated intensity of the '^CO molecule from the GRS survey data. The contours start from 15 K km s ' and the step 
is 15 K km s The velocity region included in this map is 35.1-45.7 km s The third panel shows 8 j um optical depth map o f 
the complex, constructed from the opacity maps of ~ 20 individual IRDCs identified in this region bv iPeretto & Fulleil (l2009h . 
The color scale is linear, starts from Ay k 10 mag, and the m aximu m is Ay » 90 mag. The fourth panel shows the 870 fxm dust 
emission from the region mapped by the ATLASGAL survey (ISchuller et al.ll2009l) . The color scale is linear and starts at Ay w 3.6 
mag, and the maximum is Ay « 43 mag. Lower Left: Spitzer/GLIMPSE/MIPSGAL 3-color image of the complex, showing 4.5, 
8, and 24 /im data with blue, green, and red, respectively. Lower Right: Position-velocity diagram of the region, integrated over 
-00°45' <b <Q° . The velocity interval chosen to represent the complex is indicated with dotted white lines. 
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Fig. 4. Close-up of a few IRDCs in the 35.49-00.31 complex (see Fig. lA.2l for the map of the entire complex). Left: Near-infrared 
dust extinction map. The color scale is linear and spans Ay = 3 - 30 mag. The contour shows the detection threshold of the 8 yum 
opacity data (see the right panel). Center: Integrated intensity of the ^^CO molecule in the velocity interval 24.6 - 46.9 km s '. 
Right: 8 urn optical depth map. Th e white contour is drawn at t(8 yum) = 0.3, or Ay ^6.3 mag, representing the detection threshold 
of the data (IPeretto & Fulledl2009l) . The color scale is linear with the maximum of Ay = 90 mag. 
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Fig.A.l. IRDC complex 36.74-00.16. The panels show the NIR dust extinction, '^CO line emission, 8 ;um dust opacity, and 870 
yum dust emission maps of the complex. The contours of the Ay map start from 3 mag and the step is 4 mag. The contours of the 
CO map are in multiples of 15 K km s The velocity range chosen for '^CO data is shown in the CO map. The 8 /vm color scale 
saturates white at Ay « 20 mag and the color scale of the 870 /vm map at Ay ~ 7 mag. The panels show a scale bar corresponding 
to the distance derived using foreground source count method (Dfg in Table[T]l. 
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Fig. A.2. Same as Fig. lA.ll but for IRDC complex 35.49-00.31. The 8 jim color scale saturates white at Ay ~ 60 mag and the color 
scale of the 870 fim map at Ay ~ 15 mag. See also Fig.|4]for a close-up of the central region of the complex. 
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Fig. A.3. Same as Fig. lA.ll but for IRDC complex 34.39-00.71. The 8 fim color scale saturates white at Ay ~ 40 mag and the color 
scale of the 870 fim map at Ay ~ 7 mag. 
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Fig, A.4. Same as Fig. lA.ll but for IRDC complex 34.24+00.14. The 8 jum color scale saturates white at Ay ~ 20 mag and the color 
scale of the 870 ;um map at Ay ~ 11 mag. 
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Fig. A.5. Same as Fig. lA.ll but for IRDC complex 32.09+00.09. The 8 jt/m color scale saturates white at Ay ~ 20 mag and the color 
scale of the 870 //m map at Ay ~ 7 mag. 
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Fig. A.6. Same as Fig. lA.ll but for IRDC complex 24.94-00.16. The 8 fim color scale saturates white at Ay ~ 40 mag and the color 
scale of the 870 fim map at Ay ~ 7 mag. 
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Fig. A.7. Same as Fig. lA.ll but for IRDC complex 18.54-00.16. The 8 /im color scale saturates white at Ay ~ 30 mag and the color 
scale of the 870 fim map at Ay ~ 7 mag. 
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Fig.A.8. Same as Fig. lA.ll but for IRDC complex 11.10-00.10. In addition to NIR extinction, only 8 /vm dust opacity and 
ATLASGAL data are available for this cloud. The 8 /im color scale saturates white at Ay ~ 25 mag and the color scale of the 
870 yum map at Ay ~ 7 mag. 
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Fig. A.9. Same as Fig.|Al] but for IRDC complex 03.80-01.00. In addition to NIR extinction, only ATLASGAL data are available 
for this cloud. The 8 /im color scale saturates white at Ay x 25 mag and the color scale of the 870 fim map at Ay x 7 mag. 



